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ALGORITHMS FOR SELECTION OF PRIMER PAIRS 



CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims priority to provisional application USSN 
60/317,311 filed September 5, 2001, entitled "Amplification of Nucleic Acids" which is 
incorporated by reference in its entirety for all purposes. 

COPYRIGHT NOTICE 

[0002] A portion of the disclosure of this patent document contains material 
which is subject to copyright protection. The copyright owner has no objection to the 
xerographic reproduction by anyone of the patent document or the patent disclosure exactly 
as it appears in the Patent and Trademark Office patent file or records, but otherwise reserves 
all copyrights whatsoever. 

BACKGROUND OF THE INVENTION 
[0003] The polymerase chain reaction (PGR) is a powerfiil method for amplifying 
nucleic acid sequences. Various disclosures involving this technique are found in U.S. Pat. 
Nos. 4,683,202; 4,683,195; 4,800,159; 4,965,188; and 5,512,462, each of which is 
incorporated herein by reference. In a simple form, PGR is an in vitro technique for the 
enzjmiatic synthesis of specific DNA sequences using two oligonucleotide primers that 
hybridize to complementary nucleic acid strands and flank a region that is to be amplified in a 
target DNA. A series of reaction steps of 1) template denaturation, 2) primer annealing, and 
3) extension of annealed primers by DNA polymerase, results in the geometric accumulation 
of a specific fi:-agment whose termini are defined by the 5' ends of the primers. As is well 
known, PGR is capable of selective enrichment of specific DNA sequences by a factor of 10^. 

[0004] PGR has been applied widely in molecular biology for sequencing, 
genome mapping and forensics. However, despite such wide-spread use, amplifying long 
stretches of DNA, particularly genomic DNA, is difficult. Many protocols for long range 
PGR exist; however, reaction conditions are usually optimized for amplifying specific target 
regions of interest. Applying the same "optimized" reaction conditions to amplify a different 
target region may not result in a detectable amplification product. 

[0005] In light of the above limitations, there is a need in the art for methods 
capable of amplifying nucleic acid sequences. The resulting methods may be used in some 
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embodiments to amplify mammalian target sequences across the genome to facilitate 
genotyping studies, and for other applications in the art of molecular biology. 

SUMMARY OF THE INVENTION 

[0006] The presently claimed invention provides methods for amplifying a DNA 
target sequence. One embodiment of the present invention provides robust methods for 
amplification of target sequences. In a first aspect of the invention, a method for designing 
primer pairs for the amplification reaction is provided. In a further aspect of the invention, 
reagents and cycling parameters for the amplification reaction are provided. 

[0007] Thus, the present invention provides a method for designing primer pairs 
for amplifying a target sequence, comprising the steps of: choosing a reference sequence; 
removing at least selected repeat regions in the reference sequence to yield removed and 
unremoved reference sequence; selecting primer sequences fi-om the xmremoved reference 
sequence according to two or more parameters including primer length and primer melting 
temperature yield a set of primers; evaluating the set of primers for extent of coverage and 
overlap of the reference sequence; and selecting a subset of primer pairs having reduced 
overlap from the set of primers. 

[0008] In addition, the present invention provides a method for amplifying a 
target sequence, comprising the steps of: mixing a reaction cocktail comprising 
deoxynucleotide triphosphates, target DNA, a divalent cation, DNA polymerase enzyme, a 
broad spectrum solvent, a zwitterionic buffer and at least one primer pair designed by the 
method above; heating the reaction cocktail at a denaturing temperature of about 90.0 ''C to 
about 96.0 for about 1.0 second to about 30.0 seconds; cooling the reaction cocktail at an 
annealing/extension temperature of about 50.0 "^C to about 68.0 ""C for about 1.0 minute to 
about 28.0 minutes; repeating the heating and cooling steps at least 10 times; and cooling the 
reaction cocktail to 4.0 ^'C in a final cooling step. 

[0009] Other and further objects, features and advantages would be apparent and 
eventually more readily imderstood by reading the following specification and by reference to 
the accompanying drawings forming a part thereof, or any examples of the presently 
preferred embodiments of the invention given for the purpose of the disclosure. 

DETAILED DESCMPTION OF THE FIGURES 
[0010] Figure 1 is a flow chart showing the primer pair selection process. 
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[0011] Figure 2 is a flow chart showing a detailed primer pair selection process 
according to one embodiment of the present invention. 

[0012] Figure 3 shows the sub-routines utilized to select the subset of primer pairs 
in the fourth step of the primer pair selection process. 

[0013] Figure 4 shows a basic amplification process. 

[0014] Figure 5 shows two photographs of ethidiimi bromide stained agarose gels 
on which ampUfied, genomic DNAs firom human chromosome 14 and chromosome 22 have 
been electrophoresed. 

[0015] Figure 6 shows photographs of ethidium bromide stained agarose gels on 
which amplified genomic DNA fi-om human, gorilla, chimp, and macaque has been 
electrophoresed. 

[0016] Figure 7 shows a system that may be used for designing primer pairs. 
)4p\ [0017] ^^gure 8 shows an exemplary sequence before and after removal of repeat 
(^^^quences (underline^^. 

DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION 
[0018] Reference now will be made in detail to various embodiments and 
particular applications of the invention. While the invention will be described in conjunction 
with the various embodiments and applications, it will be understood that such embodiments 
and applications are not intended to limit the invention. On the contrary, the invention is 
intended to cover alternatives, modifications and equivalents that may be included within the 
spirit and scope of the invention. In addition, throughout this disclosure various patents, 
patent applications, websites and publications are referenced. Unless otherwise indicated, 
each is incorporated by reference in its entirety for all purposes. 

[0019] The term "a" or "an" as used herein in the specification may mean one or 
more. As used herein in the claim(s), when used in conjunction with the word "comprising", 
the words "a" or "an" may mean one or more than one. As used herein "another" may mean 
at least a second or more. 

[0020] Robust methods for designing primers and amplifying target sequences are 
described herein. In one specific embodiment of the present invention, amplification of 
between about 3 kilobases and about 15 kilobases or more in length has been achieved. The 
methods result in excellent fidelity of amplification and product yield for target sequences in 
general. In some applications of the present invention, the methods result in a greater than 
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95% success rate for amplification of mammalian genomic sequences genome-wide when a 
reference sequence and a target sequence are from the same species. However, in addition, 
the methods of the present invention can be used to amplify long target sequences genome- 
wide in species closely-related to the species from which a reference sequence was taken. 
For example, human sequence can be used to design primers that will produce long-range 
amplification products of non-human primates with a success rate of greater than 80%. 
L Primer Design 

[0021] One aspect of the invention is methods for primer design. Figure 1 is a 
flow chart generally illustrating the primer selection process. In step 100 of primer design, a 
H- sequence of interest (target sequence or reference sequence) is selected for amplification and 
downloaded into a sequence file (original sequence file). The sequence file and the software 
for performing the analysis herein may be stored on a computer system such as shown in 
Figure 7. 

[0022] In step 200, repeat sequences, such as Alu and LINE sequences in the 
reference sequence, are "masked" or removed from the primer selection analysis. In step 
300, the non-repetitive, un-removed sequences that remain are analyzed according to at least 
two selection parameters and a set of all primer candidates that fit within the chosen 
parameters is established. Such selection parameters include, for example, melting 
temperature, likelihood of primer-dimer formation between the primers, primer length, and 
the like. Any of the primers generated by the third step may be used in the amphfication 
reactions of the present invention. 

[0023] In step 400, the set of primers generated by the third step is evaluated for 
coverage and overlap of the target sequence and a subset of primers is chosen so as to reduce 
the nimiber of primers needed to amplify the target sequence. 
A. Generation of a Primer Set 

[0024] In the first step 100, a sequence of interest (target sequence) may be 
obtained, for example, from public databases such as the Human Genome Project Working 
Draft team at the University of CaUfomia at Santa Cruz, NCBI, The Sanger Center, 
Whitehead Institute for Biomedical Research Center for Genome Research, Washington 
University Genome Sequencing Center, US DOE Joint Genome Institute, or Riken Gene 
Bank. Sequence generated de novo also may be used. 

[0025] The second step 200 may be performed by hand or by a computer software 
program such as, for example, the program available from the University of Washington 
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called "RepeatMasker", a program that recognizes sequences that are repeated in the genome 
(A. F. A. Smit and P, Green, 

www.genome.washington.edu/uwgc/analysistools/repeatmask, 

incorporated herein by reference). Essentially, RepeatMasker screens genomic sequences for 
repeat regions in DNA, referencing a database of known repetitive elements called RepBase. 
RepBase Version 5 has been employed in the methods of the present invention, as have 
earlier versions of RepBase. The RepBase database can be hcensed from the Genetic 
Information Research Institute (see www.girinst.org, incorporated herein by reference). 
Essentially, known repetitive sequences such as Single Interspersed Nuclear Elements 
(SDSTEs, such as alu and MIR sequences), Long Interspersed Nuclear Elements (LINEs such 
as LINEl and LINE2 sequences). Long Terminal Repeats (LTRs such as MaLRs, Retrov and 
MER4 sequences), Transposons, MERl and MER2 sequences are "masked" or removed by 
the RepeatMasker program by substituting each specific nucleotide of the repeated regions 
(A, T, G or C) with an "N" or "X". In addition, xprimer (alces.med.imm.edu. Virtual 
L Genome Center, incorporated herein by reference), a primer selection tool described below, 
can be used to identify simple, complex and intemal repeats from a small database of repeats. 
S Also, NCBI offers an Electronic PGR feature through its website (ncbi.nlm.nih.gov, 
O incorporated herein by reference). The Electronic PGR program removes repetitive 
sequences from a non-repetitive marker set. 

[0026] Figure 8 shows an exemplary sequence with repeat regions shown 
(underlined), then removed or "masked" by inserting **Ns". After the repeat regions are 
removed, primer pair candidates are selected from the unremoved sequence according to 
various parameters. 

[0027] The third step 300 may be performed by hand or by a computer software 
program. For example, commercially available software such as Primer 3 (www- 
genome.wi.mit.edu/cgi-bin/primer/primer3, incorporated herein by reference), xprimer 
(alces.med.umn.edu. Virtual Genome Center, incorporated herein by reference), Oligo 
(Molecular Biology Insights, Inc., Cascade, CO, incorporated herein by reference) or 
PrimerSelect (DNAStar, Inc., Madison, WI, incorporated herein by reference) may be 
employed. Those with skill in the art may be familiar with other programs that are available 
for primer selection or can develop such a program. In one embodiment, a software program 
is used that allows one to dictate various primer parameters such as primer melting 
temperatxu-e, primer length, stringency of hybridization, existence of duplexes, specificity of 
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hybridization, existence of a GC clamp, existence of hairpins, existence of sequence repeats, 
the dissociation minimum for a 3' dimer, the dissociation minimum for the 3' terminal 
stability range, the dissociation minimum for a minimum acceptable loop, percent maximum 
homology, percent consensus homology, the maximum number of acceptable sequence 
repeats, frequency threshold, or the maximum length of acceptable dimmers and the like. 
Also, in choosing primers for the third step, the length of a first primer of a primer pair may 
be fixed at a specific length, and the length of a second primer of the primer pair may be 
adjusted so that the melting temperature of the second primer pair is substantially the same as 
the melting temperature of the first primer. 

[0028] Primer3 is a computer program that suggests PGR primers for a variety of 
applications, for example, to create STSs (sequence tagged sites) for radiation hybrid 
mapping, or to amplify sequences for SNP discovery. Primer3 also can select single primers 
for sequencing reactions and can design oligonucleotide hybridization probes. In selecting 
ohgos for primers or hybridization probes, Primer3 can consider many factors, including 
oligo melting temperature, length, GC content, 3' stability, estimated secondary structure, the 
likelihood of annealing to or amplifying undesirable sequences (for example interspersed 
repeats), the likelihood of primer-dimer formation between two copies of the same primer, 
and the accuracy of the source sequence. In the design of primer pairs, Primer3 can consider 
product size and melting temperature, the hkelihood of primer-dimer formation between the 
two primers in the pair, the difference between primer melting temperatures, and primer 
location relative to particular regions of interest or regions to be avoided. 

[0029] xprimer is another tool for selection of PGR primers. It is designed for 
selection of sets of primers along very large queries, where the primers must all fall within a 
relatively narrow melting temperature range. It is also useful in more traditional PGR 
applications. In xprimer, the actual primer sequences are printed to standard output with some 
statistical information. At the bottom of the display, a trace shows the log probability of the 
3' end of the sequence occurring in genomic DNA as determined using a preformed database. 

[0030] PrimerSelect is a suite of tools for the design and analysis of 
oligonucleotides, including primers for PGR, sequencing, probe hybridization and 
transcription. Using DNA, RNA or back-translated proteins as templates, PrimerSelect 
details thermodynamic properties for annealing reactions. The software lists all possible 
primers, ranked in order of suitability. PrimerSelect includes a virtual lab where one can 
predict the effects the selected primers on reading fi-ames, restriction sites and other features. 
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Additionally, PrimerSelect allows for loading sequences directly from NCBI's databases, so 
that primers may be designed for published sequence. 

[0031] Oligo is a multi-functional program that searches for and selects 
ohgonucleotides from a sequence file for PGR sequencing, site-directed mutagenesis, and 
various hybridization applications. Oligo calculates hybridization temperature and secondary 
structure of oligonucleotides based on the nearest neighbor change in free energy values. 

B. Selection of a Subset of Primer Pairs 

[0032] The fourth step of primer design involves evaluating the set of primer pairs 
generated in steps one through three for coverage and overlap of the target sequence, and 
selecting a subset of primer pairs from the set of primer pairs. This fourth step may be 
performed by hand or by a computer software program. Typically the goal of the fourth step 
is to choose the primer pairs that allow one to amplify all or substantially all of the entire 
target sequence with reduced sequence amplification overlap and/or a minimal or 
substantially minimal number of primer pairs. 

[0033] In preferred embodiments, the algorithm is used to select primers that will 
ampUfy more than 90% of the unremoved target sequence, preferably more than 95% percent 
of the unremoved target sequence, and preferably more than 99% percent. Preferably the 
amplified portions of the unremoved target sequence overlap by less than 5%, preferably less 
than 2% and preferably less than 1%. Preferably a minimum or near minimum number of 
probe pairs is used. 

[0034] Algorithms known in the art may be applied for this purpose. For 
example, shortest path algorithms may be used (see, generally. Introduction to Algorithms, 
Gormen, Leiserson, and Rivest, MIT Press, 1994, pp. 514-578, incorporated herein by 
reference). In a shortest-paths problem, a weighted, directed graph G^(V,E), with weight 
function w : ^-^R mapping edges to real-valued weights is given. The weight of path p = 
(vo» vy, ....v^ is the simi of the weights of its constituent edges: 

k 

i=l 

The shortest-path weight from w to v is defined by S(u,v) being equal to min w(p):u—^v if 
there is a path from u to v, otherwise, S(u,v ) is equal to infinity. A shortest path from vertex 
u to vertex v is then defined as any path p with weight w(p) = d(u,v). Edge weights can be 
interpreted as various metrics; for example, distance, time, cost, penalties, loss, or any other 
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quantity that accumulates linearly along a path that one wishes to minimize. In the 
embodiment of the shortest path algorithm used in apphcations of this invention, each primer 
pair was considered a *Vertex". Each primer pair vertex has a relationship to each other 
primer pair vertex. This relationship is an "edge'* defined for each pair of vertices, with a 
"cost" for each edge. Cost is determined by parameters of choice, such as the extent of 
overlap of the vertices, the extent of gap between the vertices and a cost of adding another set 
of vertices to the final solution. 

[0035] Single-source shortest-paths problems focus on a given graph G=(V,E), 
where a shortest path from a given source vertex ^ e F to every vertex v e F is determined. 
Additionally, variants of the single source algorithm may be applied. For example, one may 
apply a single-destination shortest-paths solution where a shortest path to a given destination 
vertex t firom every vertex v is found. Reversing the direction of each edge in the graph 
reduces this problem to a single-source problem. Altematively, one may apply a single-pair 
shortest-path problem where the shortest path fi-om m to v for given vertices u and v is foimd. 
If the single-source problem with source vertex u is solved, the single-source shortest path 
problem is solved as well. Also, the all-pairs shortest-paths approach may be employed. In 
this case, a shortest path from w to v for every pair of vertices u and v is found — essentially, a 
single-source algorithm is run fi"om each vertex. 

[0036] One single-source shortest-path algorithm that may be employed in the 
methods of the present invention is Dijkstra's algorithm. Dijkstra's algorithm solves the 
single-source shortest-paths problem on a weighted, directed graph G=(KE) for the case in 
which all edge weights are nonnegative. Dijkstra's algorithm maintains a set of vertices, S, 
whose final shortest-path weights fi-om a source s have already been determined. That is, for 
all vertices v being elements of 5, d[v]=5(5,v). The algorithm repeatedly selects the vertex u 
as an element of V-S with the minimimi shortest-path estimate, inserts u into S, and relaxes all 
edges radiating fi^om m. In one implementation, a priority queue Q that contains all the 
vertices in V-S, keyed by their d values, is maintained. This implementation assumes that 
graph G is represented by adjacency lists. 
Dijkstra (G, w, s) 

1 INITIALIZE-SINGLE SOURCE (G,5) 

2 S^0 

3 Q^V[G] 

4 while 09^0 
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5 do EXTRACT-MiN (Q) 

6 S^S\J{u} 

7 for each vertex v e AdJ[u] 

8 do Relax (u,v,w) 
Thus, G in this case is the graph of hnear coverage of the target sequence, Q is the queue of 
all vertices to be evaluated and S is the set of vertices selected. Once one set of vertices (pair 
of primer pairs) is selected that covers a particular area of the target sequence, the other 
vertices that include these pairs can be discarded. 

[0037] Other algorithms that may be used for selecting the subset of primers 
include a greedy algorithm (again, see, Introduction to Algorithms, Gormen, Leiserson, and 
Rivest, MIT Press, 1994, pp. 329-355). A greedy algorithm obtains an optimal solution to a 
problem by making a sequence of choices. For each decision point in the algorithm, the 
choice that seems best at the moment is chosen. This heuristic strategy does not always 
produce an optimal solution. Greedy algorithms differ from dynamic programming in that in 
dynamic programming, a choice is made at each step, but the choice may depend on the 
solutions to subproblems. In a greedy algorithm, whatever choice seems best at the moment 
is chosen and then subproblems arising after the choice is made are solved. Thus, the choice 
made by a greedy algorithm may depend on the choices made thus far, but cannot depend on 
any future choices or on the solutions to subproblems. In this case, the algorithm is "greedy: 
in selecting the "best" primer pair at a moment in time according to selected criteria , without 
regard to how this selection will affect what primer pairs are available for future selection. 

[0038] One variation of greedy algorithms is Huffman codes. A Huffman greedy 
algorithm constructs an optimal prefix code and the algorithm builds a tree T corresponding 
to the optimal code in a bottom-up manner. It begins with a set of |C | leaves and performs a 
sequence of |C |-1 "merging" operations to create the final tree. For example, assuming C is a 
set of n characters and that each character c e C is an object with a defined frequency /[c], a 
priority queue keyed onf, is used to identify the two least-frequent objects to merge 
together. The result of the merger of two objects is a new object whose frequency is the sum 
of the frequencies of the two objects that were merged. For example: 

1 wHC| 

2 Q^C 

3 forz*-lton-l 

4 do z*-Allocate-NodeO 

5 X^left[z] <— EXTRACT-MIN(0 
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6 y^right[z] ^EXTRACT-MlN(0 

8 Insert (0,z) 

9 return Extract-Min(0 

Line 2 initializes the priority queue Q with the characters in C. The for loop in lines 3-8 
repeatedly extracts the two nodes x and y of lowest frequency from the queue, and replaces 
them in the queue with a new node z representing their merger. The frequency of z is 
computed as the sum of the frequencies of x and y in line 7. The node z has x as its left child 
and y as its right child. After n-l mergers, the one node left in the queue — ^the root of the 
code tree — is returned in line 9. 

P [0039] Thus, one aspect of the present invention provides a method for designing 

£ primer pairs for amplifying a target sequence, comprising the steps of choosing a reference 
% sequence; removing selected repeat regions in the reference sequence to yield removed and 
Q unremoved reference sequences; selecting primer sequences from the unremoved reference 
I' sequences according to one or more parameters to yield a set of primers; evaluating the set of 
y primers for extent of overlap and coverage of the reference sequence; and selecting a subset 
O of primer pairs having reduced overlap from the set of primers. In one embodiment of this 
S aspect of the invention, the removing step is performed by a computer program that 
references a database of known repeat sequences. Li a specific embodiment of this aspect of 
the invention, the database is RepBase. Also in a specific embodiment of the present 
invention, the computer program that performs the removing step is RepeatMasker. Another 
embodiment of this aspect of the present invention provides that one of the one or more 
parameters from the first selecting step be, for example, parameters available for selection in 
commercially-available primer selection programs such as Oligo, xprimer, PrimerSelect, 
Primer 3 and the like. Such parameters include primer melting temperature, primer length, 
stringency, existence of duplexes, specificity, GC clamp, existence of hairpins, existence of 
sequence repeats, dissociation minimum for 3' dimer, dissociation minimum 3' terminal 
stability range, dissociation minimum for minimum acceptable loop, percent maximum 
homology, percent consensus homology, maximum nxmiber of acceptable sequence repeats, 
frequency threshold, or maximum length of acceptable dimers. 

[0040] Also, in an embodiment of the present invention, the second selecting step 
selects a subset of primer pairs where this subset has a reduced number of primer pairs 
required to amplify the target sequence. Preferably, the subset is a substantially minimal 



mi 



25111369.1 
(10107967/lOllN-l) 



10 



U.S. EXPRESS MAIL#EU186312650US 



number of primer pairs required to amplify the target sequence. In one embodiment, the 
second selecting step selects the subset of primer pairs according to additional parameters 
such as length of the overlap of the target sequence amplified by the primer pairs, existence 
of gaps of target sequence between primer pairs, and the necessity of adding another primer 
pair to the subset. In an embodiment of this aspect of the invention, the second selecting step 
is performed by a computer program. Such a program may apply a shortest-paths algorithm 
or greedy algorithm, and in one embodiment of the present invention, the computer program 
applies Dijkstra's single-source shortest paths algorithm (see Figures 2 and 3). 

[0041] Figure 2 shows one embodiment of the process in Figure 1 in greater 
detail. At step 100, the target or reference sequence is downloaded from, for example, a 
public database, and stored in an original sequence file (105). At step 200, repeat sequences 
in the target sequence are removed from the primer selection process by, for example, a 
computer program such as RepeatMasker. A file of the unremoved sequence (205) is stored 
on a server or similar memory device. At step 300, primer pair candidates are selected in 
accordance with established, selected parameters, and these primer pair candidates are stored 
in a file (305) on a server or similar memory device. Preferably, all possible primer pairs that 
fall within the established parameters are stored in file 305. At step 310, the file of all 
possible primer pairs is parsed, loaded and a candidate primer pair table (315) is generated. 
At step 400, a subset of primer pairs is selected by applying, for example, a greedy algorithm. 
The subset of primer pairs is stored in file 430, a "primers to add" table, on a server or similar 
memory device. The primers to add table is then appended to a master database in step 435, 
adding this subset of primer pairs to an aggregate primer pair table 440. 

[0042] Figure 3 shows greater detail of one embodiment of step 400, selecting a 
subset of primer pairs from the table of all primer pairs generated at step 300. Step 405 
evaluates the table of all primer pairs generated at step 300, finding stretches of the target 
sequence where there are no primer pairs useful for amplification. Step 410 then adds fake 
primer pairs to cover these stretches so as to remove these gaps between primer pairs from 
the solution reached when applying the greedy algorithm in steps 415, 420 and 425. Step 415 
determines the cost of each "edge" according to pre-selected criteria for cost, step 420 finds 
the lowest cost for each set of primer pairs and step 425 finds the best path for amplifying the 
target sequence. The subset of primers generated by steps 405, 410, 415, 420, and 425 is then 
stored in a file 430 on a server or similar memory device. 
II. Computer System 
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[0043] One embodiment of the present invention provides a computer program for 
designing primer pairs for amplifying a target nucleic acid sequence. The computer program 
comprises computer code that receives input of a reference sequence; computer code that 
removes selected repeat regions in the reference sequence; computer code that selects primer 
sequences from the unremoved reference sequence; computer code that evaluates the set of 
primers for extent of coverage and overlap of the reference sequence; and computer code that 
selects a subset of primer pairs having reduced overlap from the set of primers. Preferably, 
the computer code that selects primer sequences from the unremoved reference sequence 
selects sequences according to two or more parameters including primer length and primer 
melting temperature to yield a set of primers. 

[00441 Another embodiment of the present invention provides a system that 
designs primer pairs for amplifying a target nucleic acid sequence. This system comprises a 
processor; and a computer readable medium coupled to the processor for storing a computer 
program. The computer program comprises computer code that receives input of a reference 
sequence; computer code that removes selected repeat regions in the reference sequence; 
computer code that selects primer sequences from the unremoved reference sequence; 
computer code that evaluates the set of primers for extent of coverage and overlap of the 
reference sequence; and computer code that selects a subset of primer pairs having reduced 
overlap from the set of primers. Preferably, the computer code that selects primer sequences 
from the unremoved reference sequence selects sequences according to two or more 
parameters including primer length and primer melting temperature to yield a set of primers. 

[0045] For a description of basic computer systems and computer networks, see, 
e.g.. Introduction to Computing Systems: From Bits and Gates to C and Beyond by Yale N, 
Patt, Sanjay J. Patel, 1st edition (January 15, 2000) McGraw Hill Text; ISBN: 0072376902; 
and Introduction to Client/Server Systems: A Practical Guide for Systems Professionals by 
Paul E. Renaud, 2nd edition (June 1996), John Wiley & Sons; ISBN: 0471133337, both are 
incorporated herein by reference in their entireties for all purposes. 

[0046] Appendix 1 attached hereto provides an exemplary computer code in 
Visual Basic. This code covers taking the original sequence file (105) obtained from a public 
database, through adding the subset of selected primers to the primers-to-add table (step 430) 
(see Figures 1 and 2). Figure 7 illustrates an example of a computer system that may be used 
to execute the software of an embodiment of the invention. Figure 7 shows a computer 
system 701 that includes a display 703, screen 705, cabinet 707, keyboard 709, and mouse 
711. Mouse 711 may have one or more buttons for interacting with a graphic user interface. 
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Cabinet 707 houses a floppy drive 712, CD-ROM or DVD-ROM drive 702, system memory 
and a hard drive (713) which may be utilized to store and retrieve software programs 
incorporating computer code that implements the invention, data for use with the invention 
and the Uke. Although a CD 714 is shown as an exemplary computer readable medium, other 
computer readable storage media including floppy disk, tape, flash memory, system memory, 
and hard drive may be utiUzed. Additionally, a data signal embodied in a carrier wave (e.g., 
in a network including the Internet) may be the computer readable storage medium. 
III. Amplification Reaction 

[0047] In another aspect of the present invention, methods for long range nucleic 
acid amplification are provided, including cycling temperatures, cycling times, reagents and 
reagent concentrations. The methods allow for consistent long range amplification of 
sequences genome-wide. In some embodiments of the present invention, amplification of 
between about 3 kilobases and about 15 kilobases or more in length has been achieved. In 
some applications of the present invention, the methods result in a greater than 95% success 
rate for long range amplification of mammalian genomic sequences genome-wide when the 
reference sequence and the target sequence are from the same species. However, in addition, 
the methods of the present invention can be used to amphfy long target sequences genome- 
wide in species closely-related to the species from which a reference sequence was taken. 
Various aspects of the present invention may be presented in a range format. It should be 
understood that the description in range format is merely for convenience and brevity and 
should not be construed as an inflexible limitation on the scope of the invention. 
Accordingly, the description of a range should be considered to have specifically disclosed all 
the possible subranges as well as individual numerical values within that range. For example, 
description of a range such as from 1 to 6 should be considered to have specifically disclosed 
subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 
etc., as well as individual niraibers within that range, for example, 1, 2, 3, 4, 5, and 6. This 
applies regardless of the breadth of the range. 

[0048] Figure 4 illustrates the basic steps of an amplification reaction. In step 500 
of the amplification method, reagents, target and the selected primers are combined to form a 
reaction mixture. The reaction mixture is then heated in step 505 to a temperature sufficient 
to denature the target nucleic acid, then cooled in step 510 to a temperature sufficient to allow 
annealing of the primers to the target and extension of the annealed primers. The heating step 
505 and cooling step 510 then are repeated so as to amplify the target nucleic acid. 
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[0049] Also in certain embodiments of the present invention, an initial heating 
step may be added before the heating (505)/cooling (510) cycling where the reaction cocktail 
is heated at about 90 ''C to about 96''C for 1.0 to 10.0 minutes. In a preferred embodiment, 
this initial heating step is at about 92 for about 3.0 minutes. In an alternative embodiment 
of the present invention, the cooling time for cooling step 510 may be increased for each 
successive heating/cooling cycle. In one such embodiment, the cooling time is increased by 
about 1 to about 30 seconds in each successive cycle, and in a preferred embodiment, the 
cooling time is increased by about 20 seconds in each successive cycle. 

[0050] In yet another embodiment of the present invention, an additional cooling 
step is performed after the heating (505)/cooling (510) cycle and before a final 4.0 ^^C cooling 
hold step, wherein the additional cooling step annealing/extension temperature is about 58 ""C 
to about 65 °C and is performed for about 5 minutes to about 45 minutes. In a preferred 
embodiment the additional cooling step annealing/extension temperature is about 62 °C and 
performed for about 30 minutes. 

[0051] In a specific aspect of the invention, the primers have a length of about 28 
nucleotides to about 36 nucleotides and a melting temperature of about 72.0 '*C to about 88.0 
""C. In this aspect, Tm was measured at a monovalent ion concentration of lOOOmM, a free 
Mg^ concentration of O.OmM, a total Na^ equivalent of lOOOmM, a nucleic acid 
concentration of 100 pM and where the temperature for AG calculations was 25 ""C. 

[0052] In one embodiment of the present invention, the reaction cocktail resulting 
from step 500 comprises deoxytrinucleotide triphosphates such as dATP, dTTP, dCTP, dUTP 
and dGTP or mimetics thereof, target DNA, a divalent cation, DNA polymerase enzyme, a 
broad spectrum solvent, a zwitterionic buffer and at least one primer pair designed by the 
primer selection methods described above. The heating step 505 is conducted at a denaturing 
temperature of about 90 to about 96 "^C, preferably of about 92 to about 95 ^C, and 
more preferably of about 94 °C. The denaturing temperature of the heating step 505 is 
maintained for about 1 to about 30 seconds, preferably for about 1.5 to about 5 seconds, and 
more preferably for about 2 seconds. The cooling step 510 is conducted at an 
annealing/extension temperature of about 50 ''C to about 68 **C, preferably of about 58 **C to 
about 65 ®C, and more preferably of about 62 **C. The annealing/extension temperature is 
maintained for about 1 minute to about 28 minutes, and preferably for about 15 minutes. The 
heating and cooling steps are repeated at least about 10 times and preferably about 25 to 45 
times, or more preferably about 30 to 40 times. A final cooling of the reaction cocktail to 4 
is performed after the final cooling step 510. 
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[00531 ^ embodiment of the present invention, the reaction cocktail comprises 
about 50 |iM to about 400 jxM of each primer in the primer pair, preferably about 100 nM to 
about 240 nM of each primer in the primer pair, and more preferably about 190 nM of each 
primer in the primer pair. In addition, the reaction cocktail comprises about 200 |i,M to 
about SOOjiM each dNTP, preferably about 300 ^iM to about 400p.M each dNTP, and more 
preferably about 385 [iM each dNTP. The reaction cocktail also comprises about 0.02 ng/|xl 
to about 2.0 ng/|Lil template (target) DNA, preferably about 0.05 ng/|il to about 1.0 ng/|il 
template (target) DNA, and more preferably about 0.2 ng/^1 template (target) DNA. The 
reaction cocktail may also comprise 0.0 % to about 7.0 % broad spectrum solvent, preferably 
1.5 % to about 4.5 % broad spectrum solvent, and more preferably about 3.7 % broad based 
solvent. In preferred embodiments, the broad based solvent is DMSO. 

[0054] Further, the reaction cocktail comprises 0.0 M to about 0.75 M betaine, 
preferably about 0.25 M to about 0.6 M betaine, and more preferably about 0.25 M betaine, 
and about 7 mM to about 35 mM NH4SO4, preferably about 10 mM to about 20 mM 
NH4SO4, and more preferably about 13 mM NH4SO4. The reaction cocktail also includes 
about 25 mM Tris to about 125 mM Tris, preferably about 40 mM Tris to about 80 mM Tris, 
and more preferably about 48 mM Tris, and about 100 ^iM to about 500 \iM MgCb, 
preferably about 250 ^iM to about 400 pM MgCb, and more preferably about 385 jxM MgCb. 

[0055] The reaction cocktail also comprises a polymerase. In certain 
embodiments, the reaction cocktail comprises about 0.01 units/^il to about 0.2 units/|il 
polymerase, preferably about 0.025 units/|il to about 0.07 units/^1 polymerase, and more 
preferably about 0.05 units/|al polymerase. In addition, the reaction cocktail may comprise 
about 0 mM to about 50 mM zw^itterionic buffer, preferably about 10 mM to about 30 mM 
zwitterionic buffer, and more preferably about 25 mM zwitterionic buffer. In some 
embodiments, the zwitterionic buffer is Tricine. 

[0056] Also in some embodiments, about 0.005 iig/\xl to about 0.10 ^g/[il taq 
antibody may be added to the reaction cocktail. Preferably, about 0.01 ^g/|iil to about 0.05 
ixg/^il taq antibody is added to the reaction cocktail, and more preferably about 0.025 ^g/|xl 
taq antibody is added to the reaction cocktail. 
IV, Applicability to Diverse Sequences 

[0057] PGR has been applied widely in molecular biology; however, despite such 
wide-spread use, amplifying varying long stretches of DNA is difficult. Many protocols for 
long range PGR exist; however, reaction conditions are usually optimized for amplifying 
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specific target regions of interest. Similar amplification success is not achieved when these 
"optimized" reaction conditions are used on different target regions. In the present invention, 
however, amplification of between about 3 kilobases and about 15 kilobases or more in 
length has been achieved on varied genomic sequences genome-wide. The methods result in 
excellent fideUty of amplification and product jdeld for mammalian target sequences in 
general. In some applications of the present invention, the methods result in a greater than 
95% success rate for amplification of mammalian genomic sequences when the reference 
sequence and the target sequence are firom the same species. However, in addition, the 
methods of the present invention can be used to amplify long target sequences genome-wide 
in species closely-related to the species from which a reference sequence was taken. For 
p example, human sequence can be used to design primers that will produce long-range 
iS amplification products of non-human primates with a success rate of greater than 80%. 

[0058] Figure 4 shows the results obtained with the methods of the present 
invention for human chromosome 14 sequence used as a reference sequence for primer 
design and human target DNA and human chromosome 22 sequence used as a reference 
sequence for primer design and human target DNA. Figure 5 shows the results obtained with 
the methods of the present invention with human DNA used as a reference sequence for 
^5 primer design and human, gorilla, chimpanzee, and macaque genomic DNA used as target 
sequences. 
V. Examples 

[0059] The examples below illustrate specific implementations of the inventions 
described herein. 

A. Preparation and Scoring of Somatic Cell Hybrids 

[0060] Standard procedures in somatic cell genetics were used to separate human 
DNA strands (chromosomes) fi-om a diploid state to a haploid state. Diploid hxmian 
lymphoblast cell lines from a human diversity panel lymphoblast line (available fi-om Coriell 
Cell Repositories, Camden, NJ) were fused to a diploid hamster fibroblast cell line containing 
a mutation in the thymidine kinase gene. In a sub-population of the resulting fused cells, 
human chromosomes were introduced into the hamster calls. Selection for the human DNA- 
containing hamster cells (fusion cells) was achieved by utilizing HAT medium. Only 
hamster cells that had a stably incorporated human DNA strand grow in cell culture medium 
containing HAT. 



fy 
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[0061] Hamster cell line A23 cells were pipetted into a centrifuge tube containing 
10 ml DMEM in which 10% FBCS + IX Pen/Strep + 10% glutamine were added, 
centrifuged at 1500 rpm for 5 minutes, resuspended in 5 ml of RPMI and pipetted into a 
tissue culture flask containing 15 ml RPMI medium. The lymphoblast cells were grown at 
3T C to confluence. At the same time, human lymphoblast cells were pipetted into a 
centrifuge tube contaming 10 ml RPMI in which 15% FBCS + Ix Pen/Strep -h 10% glutamine 
were added, centrifuged at 1500 rpm for 5 minutes, resuspended in 5 ml of RPMI and 
pipetted into a tissue culture flask containing 15 ml RPMI. The lymphoblast cells were 
grown at ST'^C to confluence. 

[0062] To prepare the A23 hamster cells, the media was aspirated and the cells 
were rinsed with 10 ml PBS. The cells were then trypsinized with 2 ml of trypsin and 
divided into 3-5 plates of fresh media (DMEM without HAT) and incubated at The 
lymphoblast cells were prepared by transferring the culture into a centrifuge tube and 
centrifuging at 1500 rpm for 5 minutes, resuspending the cells in 5 ml RPMI and pipetting 1 
to 3 ml of cells into 2 flasks containing 20 ml RPMI. 

[0063] To achieve cell fusion, approximately 8-10 x 10^ lymphoblast cells were 
centrifuged at 1500 rpm for 5 min. The cell pellet was then rinsed with DMEM by 
resuspending the cells and centrifuging them again. The lymphoblast cells were then 
resuspended in 5 ml DMEM. The recipient A23 hamster cells had been grown to confluence 
and split 3-4 days before the fusion and were, at this point, 50-80% confluent. The old media 
was removed and the cells were rinsed 3 times with DMEM and finally suspended in 5 ml 
DMEM. The lymphoblast cells were slowly pipetted over the recipient A23 cells and the 
combined culture was swirled slowly before incubating at 2>TC for 1 hour. After incubation, 
the media was gently aspirated from the A23 cells, and 2 ml room temperature PEG 1500 
was added by touching the edge of the plate with a pipette and slowly adding PEG to the 
plate while rotating the plate with the other hand. It took approximately one minute to add all 
of the PEG in one full rotation of the plate. Next, 8 ml DMEM was added down the edge of 
the plate while rotating the plate slowly. The PEB/DMEM mixture was aspirated gently from 
the cells and then 8 ml DMEM was used to rinse the cells. This DMEM was removed and 10 
ml fresh DMEM was added and the cells were incubated for 30 min. at 37**C. Again the 
DMEM was aspirated from the cells and 10ml DMEM in which 10% FBCS and Ix Pen/Strep 
were added, was added to the cells, which were then allowed to incubate overnight. 

[0064] After incubation, the media was aspirated and the cells were rinsed with 
PBS. The cells were then trypsinized and divided among 20 plates containing selection 



25111369.1 
(10107967/101 lN-1) 



U.S. EXPRESS MAIL#EU186312650US 



media (DMEM in which 10% FBCS + Ix Pen/Strep + Ix HAT were added) so that each plate 
received approximately 100,000 cells. The media was changed on the third day following 
plating. Colonies were picked and placed into 24-well plates upon becoming visible to the 
naked eye (day 9-14). If a picked colony was confluent within 5 days, it was deemed healthy 
and the cells were trypsinized and moved to a 6-well plate. 

[0065] DNA and stock hybrid cell cultvires were prepared from the cells from the 
6-well plate cultures. The cells were trypsinized and divided between a 100 mm plate 
containing 10 ml selection media and an eppendorf tube. The cells in the tube were pelleted, 
resuspended 200 ^il PBX and DNA was isolated using a Qiagen DNA mini kit at a 
concentration of <5 million cells per spin column. The 100 mm plate was grown to 
confluence, and the cells were either continued in culture or frozen, 

[0066] Scoring for the presence, absence and diploid/haploid state of each hybrid 
was performed using the Affymetrix, Inc. HuSNP GENECHIP® (Affymetrix, Inc. of Santa 
Clara, CA, GENECHIP® HuSNP Mapping Assay, reagent kit and user manual, Affymetrix 
Part No. 900194), which can score 1494 markers in a single chip hybridization. As a control, 
the human diploid lymphoblast cell line was screened using the HuSNP chip hybridization 
assay, and any SNPs which were heterozygous in the parent lymphoblast diploid cell line 
were scored for haploidy in each fusion cell line. By comparing the markers that were 
present as "AB" heterozygous in the parent diploid cell line to the same markers present as 
"A" or "B" (hemizygous) in the hybrids, the human DNA strands which were in the haploid 
state in each hybrid line was determined. 

B. Primer Selection 

[0067] Human genomic sequence was used as a reference sequence for primer 
selection in this example of the present invention, and human genomic DNA derived from 
somatic cell hybrids was used as target DNA. In addition, in an altemative appUcation of the 
present invention, himian genomic sequence was used as reference sequence for primer 
selection and genomic DNA from gorilla and chimpanzee was used as target DNA. 

[0068] Figure 2 is a flow chart showing a detailed primer selection process 
according to one embodiment of the present invention. The first step 100 of primer selection 
required selecting a sequence of interest (target sequence or reference sequence) and creating 
an original sequence file (105) containing this selected sequence. Next, repeat regions in the 
target sequence were removed (200), and a removed file was created containing the 
unremoved sequence (205). In the third step, the sequences in the removed file were run 
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through a primer pair selection program (300), and the set of all possible primers generated 
was stored in an oligo output file (305). The information from the oligo output file was then 
used to create a candidate primer pair table (315). In step four of the selection process (400), 
an optimal subset of primer pairs was selected from the set of all possible primer pairs in the 
primer pair table. The output from the selection of the optimal subset of primer pairs was 
stored in the primers to add table (430), which was then appended to the master database 
(435) and stored in an aggregate primer pair table (440). 

[0069] First, human sequence to be used as the reference sequence for primer 
design was acquired from the Human Genome Project Working Draft team from the 
University of Califomia at Santa Cruz where sequence assembly was performed using 
sequences obtained from the High Throughput Genomic Sequence (HTGS) database. The 
HTGS database is a public database with sequences contributed by, inter alia, the Human 
Genome Project Working Draft team. The UTSC assembly is available at the UCSC site 
[http://genome.cse.ucsc.edu/], and a detailed description of the data format can be found at 
[http://genome.cse.ucsc.edu/goldenPath/datorg,html]. Sequence was also acquired from 
NCBI. 

[0070] In the second step, acquired reference sequence was processed by a 
software program called "RepeatMasker", available for licensing from the University of 
Washington (see:A. F. A. Smit and P. Green, 

[www.genome.washington.edu/uwgc/analysistools/repeatmask.htm]). 

RepeatMasker screens genomic sequences for repeat regions in DNA, referencing a database 
of known repetitive elements called RepBase. RepBase Version 5 was employed in the 
methods of the present invention, as were earlier versions of RepBase, The RepBase 
database was licensed from the Genetic Information Research Institute (see www.girinst.ofg). 
Known repetitive sequences such as Single Interspersed Nuclear Elements (SINEs, such as 
alu and MIR sequences). Long Interspersed Nuclear Elements (LINEs such as LINEl and 
LINE2 sequences). Long Terminal Repeats (LTRs such as MaLRs, Retrov and MER4 
sequences), Transposons, MERl and MER2 sequences were "masked" or removed by the 
RepeatMasker program by substituting each specific nucleotide of the repeated regions (A, T, 
G or C) with an "N" or "X". Local nucleotide duplications were not removed. In one 
application of the present invention, the default settings of RepeatMasker were used, and the 
human.ref \\hr?ccy (human repetitive elements) and simple.ref were concatenated and 

combined to SnRNAs from the pseudo.refWbxdiry to create a "custom" library. Those skilled 
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in the art will appreciate that any computer program, algorithm or selection process, 
including manual selection, which identifies and eliminates from primer selection repetitive 
sequences from the reference sequence may be used as an alternative to RepeatMasker. 

[0071] Once the reference sequence was masked and repetitive regions removed, 
a third step was performed where the unremoved sequence output was then entered into the 
commercially-available primer design program, Oligo 6.52 using the following search 
parameters: 

Search for: Primers and Probes 

j^Strand Search 
Select: 

Complex Substrate 

Compatible Pairs 

Duplex-free Oligonucleotides 

Highly Specific Oligos [3 '-end stability] 

Oligonucleotide with GC Clamp 

Eliminate False Priming Oligonucleotides 

Oligonucleotides within Selected Stability Limits 

Hairpin-free Oligonucleotides 

Eliminate Homooligomers/Sequence Repeats 

Eliminate Frequent Oligos 

Search Mode: Mark 

PCR Product Length: 3000 to 15000 

General Settings: 

High Search Stringency 

No Auto Change 

Adjust Length to Match Tm's 

Parameters: 

Oligonucleotide Length: 32 nt 
Acceptable 3'-Dimer AG: -3.5 kcal/mol 
Maximum Length of Acceptable Dimers: 4 Base Pairs 
3 '-terminal Nucleotides Checked for Dimers: 23 
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3'-terminal Stability Range: -5.5 to -9.8 kcal/mol 

GC Clamp Stability: -10.0 kcal.mol 

Minimum Acceptable Loop AG: 0.0 kcal/mol 

Oligo Tm Range [58.1 to 108.1]: 72.0 to 88.0 ^'C 

Max Acceptable False Priming Efficiency: 170 Points 

Min Consensus Priming Efficiency: 340 Points 

Max Acceptable Homology: 50% 

Min Consensus Homology: 95% 

Max Number of Acceptable Sequence Repeats: 3 

Max Degeneracy: 1 

Frequency Threshold: 1000 

Non-Search Parameters: 

Monovalent Ion Concentration: lOOOmM 
Free Mg^ Concentration: O.OmM 
Total Na"^ Equivalent: lOOOmM 
Nucleic Acid Concentration: 100 pM 
Temperature for AG Calculations: 25 **C 

All possible primer pairs generated within the established parameters were saved to a file. 
Any of the generated primer pairs may be used in the amplification reactions of the present 
invention; however, typically primer pairs will be chosen that cover as much of the reference 
sequence as possible with reduced overlap. 

[0072] In the present embodiment, the primer pair set output obtained firom Oligo 
6.52 was, in the fourth step of primer selection, subjected to Dijkstra's algorithm (again, see 
Introduction to Algorithms, Cormen, Rivest and Leiserson (1990); ISBN 0262031418)), The 
goal of this step being to find a best subset of primer pairs to amplify the target sequence out 
of all possible sets of primer pairs generated by Oligo 6.52. Dijkstra's algorithm solves the 
single-soiurce shortest path problem on a weighted, directed graph. In the embodiment of this 
algorithm used in applications of this invention, each primer pair was considered a "vertex" 
with an "edge" defined for each pair of vertices. An associated "cost" was assigned to each 
edge where the cost reflected the amount of: 1) the overlap of vertices (cost = the length of 
the overlap); 2) the gap between two primer pairs (cost = lOx the length of the gap); and 3) a 
fixed value for having to add another vertex to the set (which increased the number of 
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primers that must be used) (cost for additional primer pair = 4000). In one application of the 
present invention, the path with the lowest cost was selected, where total cost equals the sum 
of the costs of edges in the path. For example, assume three exemplary primer pairs: 

5' position of the forward primer 5* position of the reverse primer 

Primer 1: 1000 2000 
Primer 2: 1800 3000 
Primer 3: 2100 4000 

The "edges" are defined as being between Primer 1 and Primer 2, Primer 1 and Primer 3, and 
Primer 2 and Primer 3. The cost associated with the edge Primerl/Primer2 is 200 + 0 (100) + 
4000 = 4200 (reflecting the 200 base overlap between the amplicons). The cost associated 
with edge Primerl/Primer3 is 0 + 10 (100) + 4000 = 5000 (reflecting the 100 base pair gap 
between Primer 1 and Primer 3). The cost associated with edge Primer2/Primer 3 would be 
900 + 0 (100) + 4000 = 4900 (reflecting the 900 base overlap between the amplicons). 

[0073] In one embodiment of the present invention, the computer code for 
evaluating the primer set for extent of coverage and overlap of the target sequence and 
selecting the subset of primer pairs was comprised of a main module, a first level subroutine, 
and several second level subroutines. This code is reproduced below. The main module. 
Main, essentially loaded the file of all possible primer pairs from the imremoved sequence 
from the third step, ran an error check on the sequences to assure the primers pairs were 
unique, ran the first level subroutine, then took the information output from the first level 
subroutine and appended this information to a local repository of information, which 
ultimately was copied to an aggregate primer pair table. 

[0074] The first level subroutine, Select Optimal Primers, directed several second 
level subroutines, which essentially applied Dijkstra's algorithm to select a subset of primer 
pairs from the set of all possible primer pairs (see Figure 3). Select Optimal Primers 
retrieved the information from the primer pair table (Parse Oligo Results Files), foimd gaps in 
the primer pair amplification coverage of the target reference sequence (Find Gaps 405), 
added fake primer pairs to cover the gaps so as not to penalize the solution for the subset 
selection for an xmavoidable gap (Add Fake Primer Pairs for Gaps 410), determined a cost for 
each edge (Find Edges 415), computed the lowest cost for every possible set of primer pairs 
(Compute Minimxmi Costs 420), found the best subset of primer pairs (Find Best Path 425), 
and added this subset of primer pairs to a local repository (430) that was then added to the 
final aggregate repository of primer pairs (440). 
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C. Amplification Reaction 

[0075] The amplification reaction involves both an amplification reaction mix or 
cocktail and thermocycling parameters. In one application of the present invention, the 
reaction mix was prepared by making two master reaction mixes, then adding an aliquot of 
each mix to the primer pairs in the following manner: 
PCR setup: 

13 nL total volume reactions 

Master Mix 1 : 

Amount per reaction Final Concentration per reaction 



Reagents: 

Water 

dNTPs, lOmM each 
template DNA (20ng/ [iL) 
10% DMSO/5M betaine 
Total Volume: 



4.775 
0.5 \\L 
0.1 ^iL 
.625 ^iL 
6 ^iL 



385 ^iM 

2ng 

0.48X 



Master Mix 2: 

Water 3.5625 [aL 

140 mM NH4SO4/5OO mM Tris 1 .25 

25 mM MgCl2 2.7 \xL 

Taq Polymerase (2.5U/ |iL) 0.2625 |iL 

DMSO 0.4 jiL 

Tricine (IM) 0.325 |aL 

Total Volume: 6.0 ^iL 



13mM/48mM 
385 |LiM 
0.66 imits 
3.1% 
25 mM 



The Master Mixes were prepared and kept on ice. 6.0 |jL of each Master Mix was added to 
tubes containing 1 |iL of the primers where the primers contained 2.5 ^M of each of the 
forward and reverse primers for a final concentration of 192 nM each primer in the final 13 
|liL reaction volimie. 

[0076] In an altemative embodiment of the present invention, the taq polymerase 
can be eliminated from Master Mix 2, and instead combined with 0.015|Lig/|iL TaqStart 
antibody and buffer to form an antibody-bound taq complex which is then added to the 
reaction cocktail after the Master Mix 1 and 2 have been combined. 
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[0077] Reagents for the reaction cocktails can be obtained from the following 
sources: dNTP's (Life Technologies), Taq polymerase (Roche Molecular Biosciences, 
Epicentre Tecnlogies, Biorad Laboratories or Applied Biosystems), tricine, tris, NH4SO4, 
MgCl2,betaine, and DMSO (Sigma Aldrich), Taqstart antibody (Clontech), 

[0078] In one example, the cycling conditions were as follows: 

Initial heating step: 94®C for 3 minutes 

10 cycles of: heating step: 94**C for 2 seconds 

cooling step: 62*'C for 15 minutes 

28 cycles of: heating step: 94®C for 2 seconds 

cooling step: 62°C for 15 minutes for the first cycle, with an increase in time 
of 20 seconds in each subsequent cycle 

B Final cooling step: 62°C for 25 minutes 

= 4»C hold 

Q Also, in an altemative example of the present invention, the cycling conditions were as 

P follows: 

jp^ Initial heating step: 94**C 3 minutes 

H 35 cycles of: heating step: 94**C for 2 seconds 

cooling step: 62®C for 12 minutes 
y Final cooling step: 62**C for 25 minutes 
' 4«C hold 

Aliquots of each completed amplification reaction were run on a 0.8% agarose gel and 
visualized with ethidium bromide. 

[0079] The above description is illustrative and not restrictive. Many variations 
of the invention will become apparent to those of skill in the art upon review of this 
disclosure. The scope of the invention should, therefore, be determined with reference to the 
appended claims along with their full scope of equivalents. 
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APPENDIX 1 

© Copyright, Perlegen Sciences, Inc., All Rights Reserved 

MAIN MODULE: 

MAIN - Parse And Select Primers - this routine drives everything 



Upon start-up, run this routine. Check the table "PPC" to find out 
whether 

to start a contig. (Populated "PPC" <==> don't start). If "PPC" is 
empty, 

• this process will attempt to process a batch of primer pairs (for a 
contig) , 

» It will do these steps: 





(A) 


Parse the OLIGO files for this contig 


- > 


table 


"PPC" (see below) 




(B) 


Select unique primer pairs from "PPC" 


-> 


table 


"PrimerPair" 




(C) 


Run "SelectOptimalPrimerPairs" : PP ' s 


- > 


table 


" PrimersToAdd" 


3 


(D) 

I 


Append primer pairs from PrimersToAdd 


- > 


table 


"Primers" (below) 



knows 
tables : 



This database is assumed to have a contig for its name. That's how it 
which contig to do. Also, it needs to have the following linked 



(1) Chrlnfo 

(2) Ctglnfo 

(3) Primers 



Information specific to this chromosome 
Lengths of all contigs for this chromosome 
Where to append the selected primer pairs 



I I I I t I I T I I I f I I I I t t f t I I I I f I 1 1 I f I I I f I I I I 1 I I I I I f t I I I I I I I I I t I I t f I I I I I I I I I I I I I 



Option Compare Database 
Option Explicit 



* Only Functions can be called in AutoExec Macro! 
Public Function ParseAndSelectPrimers ( ) As Boolean 

Dim rst As Recordset 

Dim strPath As String 

Dim strContig As String 

Dim strPrefix As String 

Dim IngContigLen As Long 

Dim IngSegmentLen As Long 

Dim IngOverlap As Long 
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Dim IngSeqlD As Long 
Dim NumPPs As Long 

' Find out whether to do anything : 
Set rst = CurrentDb.OpenRecordset ("PPC") 
If rst .RecordCount <> 0 Then 
DoCmd . OpenForm "frmNotes" 
Exit Function 
End If 

Set rst = Nothing 
' Initialize: 

strPath = GetltemValue ("OligoResultsPath" ) 
strContig = GetContig 

strPrefix = GetltemValue ( "FileNamePrefix" ) 
IngContigLen = GetContigLen ( strContig) 
IngSegmentLen = GetltemValue ( "SegmentLength" ) 
IngOverlap = GetltemValue ( "Overlap" ) 
IngSeqlD = GetltemValue ( "SeqID" ) 

' (A) 

WriteLog "STARTING..." 

WriteLog "ParseOligoFileSet " & strPath & strContig & " - Length: 
" & IngContigLen 

NumPPs = ParseOligoFileSet (strPath, strContig, strPrefix, 
IngContigLen, IngSegmentLen, IngOverlap, IngSeqlD) 

• (B) 

WriteLog NumPPs & " primer pairs found. Append unique ones to 
PrimerPair . . . " 

Cur rentDb. Execute "Append PPC -> PrimerPair" 

• (C) 

WriteLog " SelectOpt imalPrimerPairs " 
SelectOptimalPrimerPairs 

» (D) 

WriteLog "Append PTA -> Primers" 
CurrentDb. Execute "Append PTA -> Primers" 
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' That ' s it : 

WriteLog "Application Quit - " & strContig 
Application . Quit 
End Function 



FIRST LEVEL SUBROUTINE; SELECT OPTIMAL PRIMER PAIRS 

option Compare Database 
Option Explicit 



' Find optimal bunch of primer pairs : 

N- » Assume PrimerPair is ready (local) and Edge is not indexed and Gap 

M exists 

£ Public Sub SelectOptimalPrimerPairs ( ) 

Fy WriteLog "FindGaps" 

S FindGaps 

1,, WriteLog "AddFakePrimerPairsToCoverTheGaps " 

14: AddFakePrimerPairsToCoverTheGaps 

£5 WriteLog "FindEdges" 

Hi FindEdges 



WriteLog "Createlndexes" 
IndexFieldlnTable "Src", "Edge" 
IndexFieldlnTable "Dst", "Edge" 
IndexFieldlnTable "Cost", "Edge" 



WriteLog "Executing queries - ZeroOut COST, etc." 
CurrentDb . Execute "ZeroOut COST, PRED, DONE" 
CurrentDb . Execute "ZeroOut COST of PPO" 



WriteLog "ComputeMinCosts " 
Comput eMinCos t s 

WriteLog "Initialize field ^SELECTED'" 
RenameFieldlnTable "DONE", "SELECTED", "PrimerPair" 
CurrentDb. Execute "UPDATE PrimerPair SET SELECTED = No" 
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WriteLog "FindBestPath" 
FindBestPath 

WriteLog "Queries - Make Selected, PrimersToAdd" 
CurrentDb. Execute "Make Selected" 
CurrentDb . Execute "Make PrimersToAdd" 

WriteLog "Skipping - FindActualGaps - can run on main machine!" 
* Skip - FindActualGaps 

WriteLog "FINISHED!" 
y^.-. MsgBox "Primer Pair Optimizer finished!", vbinf ormation. Now 

2 End Sub 

SECOND LEVEL SUBROUTINES: 

J: 

r 1. PARSE OLIGO RESULTS 

fkJ Parse Oligo Results File(s): 

I ss = = ss = =: = = = = = = = = = = = = s = =: = =s = ssssss 

' Get primer pairs from Oligo results files. Store them locally (table 
p "PPC") . 

' Specifically, get ALL the primer pairs for a chromosome (many 

contigs) . 

' Assumptions: (1) OLIGO results filenames are like [Startbase] . txt 
I 

I r I 1 I I I t I I I I I I I I t I I I I I I I I I f f I f t I I I I t t I I I I I t I I I t I I I I I I I I I I I I I I f t I t I t I I t I I f 

Option Compare Database 
Option Explicit 

' Parse an OLIGO results file SET (all for one contig) : 

Public Function ParseOligoFileSet (ByVal strOligoParentPath As String, 
ByVal strContig As String, _ByVal strFileNamePref ix As String, ByVal 
IngContigLen As Long, _ ByVal IngSegmentLen 

As Long, ByVal IngOverlap As Long, _ 

ByVal IngSeqlD As Long) As Long 

Dim I As Long 

Dim strFileName As String 

' Initialize: 
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ParseOligoFileSet = 0 

• Parse the primer pairs : 

For I = 1 To IngContigLen - IngOverlap Step IngSegmentLen - 

IngOverlap 

strFileName = strOligoParentPath & strContig & "\" & CStr(I) 

& " . txt " 

Say "Primer pairs found: " & ParseOligoFileSet & ": 

Parsing " & strFileName 

ParseOligoFileSet = ParseOligoFileSet + 

ParseOligoFile (strFileName, strContig, I, strFileNamePref ix, IngSeqlD) 

Next I 

' Done : 
Say "Ready" 
End Function 

• Write primer pairs to table "PPC" : 



3 



Private Function ParseOligoFile (ByVal strFileName As String, ByVal 
StrContig As String, ByVal IngStartBase As Long, _ByVal StrFileNamePref ix 
rf As String, ByVal IngSeqlD As Long) As Long 

O Dim rst As Recordset 

^ Dim iFileNum As Integer 

jfij Dim strLme As String 

Dim IngPairNum As Long 

Dim IngPrimerLen As Long 

Dim nColonPosn As Long 

Dim nLetterPosn As Long 

Dim nThreePosn As Long 

' Open the table : 

Set rst = CurrentDb.OpenRecordset ( "PPC" ) 
With rst 

• Open the file: 
iFileNum = FreeFile 

Open StrFileName For Input As # iFileNum 

' Verify that the file ID matches the file name in line 2 of 
the Oligo file: 

Line Input #iFileNum, strLine 
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Line Input #iFileNum, strLine 

If IngStartBase <> CLng (GetSubstring ( 1 , strLine, 
strFileNamePref ix, "_")) Then 

MsgBox IngStartBase & " not found in " & strLine, 
vbCritical, "Possible Parsing Problem" 

Stop 

End If 



' Get all the primer pairs : 
Do Until EOF (iFileNum) 



S?— lis 



2V' 



* Input a line from the file: 
Line Input #iFileNum, strLine 

* Check for new primer pair: 

If Left$ (strLine, 6) = "Pair #" Then 

* Add the new pair: 
. AddNew 

iSequencelD = IngSeqlD 

IContig = strContig 

IFilelD = IngStartBase 

IngPairNum = CLng (iyiid$ (strLine, 7) ) 

! PairNum = IngPairNum 



' Product Length : 
Line Input #iFileNum, strLine 
lAmpliconLen = Val (MidS (strLine, 16)) 



' Forward Coordinates : 
SkipLines 3, iFileNum 
Line Input #iFileNum, strLine 
nColonPosn = InStr (strLine, ":") 
nLetterPosn = InStr (strLine, "U") 

!FPOS = IngStartBase + CLng (Mid$ (strLine, nColonPosn 
+ 1, nLetterPosn - (nColonPosn + 1) ) ) -1 

! ForwardLen = Val (Mid$ (strLine , nLetterPosn + 1)) 



' Forward Sequence : 

Line Input #iFileNum, strLine 
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nThreePosn = InStr (strLine , "3") 

1 ForwardSeq = RemoveWhiteSpace (Snip (strLine , 4, 



nThreePosn) ) 



• Forward Tm: 

Line Input #iFileNum, strLine 
IForwardTm = Val (Mid$ (strLine , 4)) 



• Reverse Coordinates : 
SkipLines 4, iFileNum 
Line Input #iFileNum, strLine 
nColonPosn = InStr (strLine , ":") 
O nLetterPosn = InStr (strLine, "L") 

M IngPrimerLen = Val (Mid$ ( strLine , nLetterPosn + 1)) 

p|H IREND = IngStartBase + CLng (Snip (strLine, nColonPosn 

j| +1, nLetterPosn)) + IngPrimerLen - 2 

l5 IReverseLen = IngPrimerLen 

w - 
3 

' Reverse Sequence : 
Line Input #iFileNum, strLine 
nThreePosn = InStr (strLine , "3") 

IReverseSeq = RemoveWhiteSpace (Snip ( strLine , 4, 



nThreePosn) ) 



* Reverse Tm: 

Line Input #iFileNum, strLine 
IReverseTm = Val (Mid$ (strLine, 4)) 



' Update: 
. Update 
End If 

Loop 
End With 

Set rst = Nothing 
Close #iFileNum 

' Return the number of pairs parsed (i.e., the most recent pair 

) : 

ParseOligoFile = IngPairNum 
End Function 
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2. FIND GAPS 




ion Compare Database 
ion Explicit 

ilic Sub FindGapsO 

Dim IngMaxREND As Long 

Dim IngREND As Long 

Dim IngFPOS As Long 

Dim sSQL As String 

Dim rstMain As Recordset 

Dim rstGap As Recordset 



Say "Finding the gaps. . . " 

sSQL = "SELECT * FROM PrimerPair order by FPOS, PPCID" 
Set rstMain = Cur rent Db . OpenRecordset (sSQL) 
Set rstGap = CurrentDb . OpenRecordset ( "Gap" ) 
With rstMain 

IngMaxREND = 0 

Do Until .EOF 



IngREND = I REND 
IngFPOS = i FPOS 

' Check for gap: 

If IngMaxREND + 1 < IngFPOS Then 
With rstGap 



. AddNew 

ILastREND = IngMaxREND 
!NextFPOS = IngFPOS 
iNextPPCID = rstMain I PPCID 
- Update 



End With 
End If 

' Update Max REND: 

If IngMaxREND < IngREND Then 

IngMaxREND = IngREND 
End If 
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• status : 

If . AbsolutePosition Mod 1000 = 0 Then 

Say "Finding Gaps: " & .AbsolutePosition 
End If 

. MoveNext 

Loop 

' Fake a gap from Max REND to ONE BILLION: 
With rstGap 
. AddNew 

! La St REND = IngMaxREND 
INextFPOS = 1000000000 
. Update 



B End With 

p End With 



Set rstMain = Nothing 
Set rstGap = Nothing 
Say "Ready" 



End S\ib 



3. ADD FAKE PRIMER PAIRS FOR GAPS 

Option Compare Database 
Option Explicit 

Public Sub AddFakePrimerPairsToCoverTheGaps ( ) 
Dim IngPPCID As Long 
Dim rstGap As Recordset 
Dim rstPP As Recordset 

Set rstGap = CurrentDb . OpenRecordset ( "SELECT * FROM Gap ORDER BY 

LastREND") 

Set rstPP = CurrentDb. OpenRecordset ("PrimerPair") 
Say "Adding fake PP's to cover the gaps..." 

IngPPCID = 0 ' Woe to you who changes this 

line 

Do Until rstGap. EOF 
With rstPP 
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already) 



. AddNew 

IPPCID = IngPPCID 
IFPOS = rstGap ILastREND 
IREND = rstGap INextFPOS 
.Update 
End With 

IngPPCID = IngPPCID - 1 ' Avoid positive ID*s (used 

rstGap . MoveNext 

Loop 

Set rstGap = Nothing 
Set rstPP = Nothing 
Say "Ready" 



End Sub 



4. FIND EDGES 

option Compare Database 
Option Explicit 

* How much worse is a skipped base than an overlapping base: 
Private Const mcnGapPenaltyPerBase As Long = 10 

' How much worse is an additional amplicon than an overlapping base: 
Private Const mcnPenaltyPerAmplicon As Long = 4000 



' Populate the Edge table with costs of amplicon pairing: 
Public Sub FindEdgesO 

Dim rstSrc As Recordset 

Dim rstDst As Recordset 

Dim rstEdge As Recordset 

Dim sSQL As String 

Dim IngSrcID As Long 

Dim IngSrcFPOS As Long 

Dim IngSrcREND As Long 

Say "Finding Edges - Initializing Recordset..." 
sSQL = "SELECT * FROM PrimerPair ORDER BY FPOS, REND" 
Set rstSrc = CurrentDb .OpenRecordset (sSQL) 
Set rstEdge = CurrentDb .OpenRecordset ( "Edge" ) 
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Do Until rstSrc.EOF 

IngSrcFPOS = rstSrclFPOS 
IngSrcREND = rstSrclREND 
IngSrcID = rstSrclPPCID 
If IngSrcID Mod 11 = 0 Then 

Say "Finding Edges for " & IngSrcID 
End If 

sSQL = "SELECT * FROM PrimerPair WHERE " & _ 
" FPOS > " & IngSrcFPOS & " AND " & _ 
" REND > " & IngSrcREND & " AND " & 
" FPOS < " & IngSrcREND + 100 0 & _ 
" ORDER BY FPOS, REND " 
Set rstDst = CurrentDb . OpenRecordset ( sSQL) 
Do Until rstDst. EOF 
With rstEdge 
O . AddNew 

" ! Src = IngSrcID 

g !Dst = rstDst !PPCID 



O 



iCOST = GetCost (IngSrcREND, rstDst! FPOS) 
. Update 



End With 
rstDst .MoveNext 

Loop 

rstSrc .MoveNext 

Loop 

Set rstSrc = Nothing 
Set rstDst = Nothing 
Set rstEdge - Nothing 
Say "Ready" 



End Sub 



5, COMPUTE MINIMUM COSTS 

Option Compare Database 
Option Explicit 

Public Sub ComputeMinCosts ( ) 

Dim IngMinSrcCost As Long 
Dim IngSrcID As Long 
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Dim IngEdgeCost As Long 
Dim sSQL As String 
Dim rstMin As Recordset 
Dim rstSrc As Recordset 
Dim rstEdge As Recordset 
Dim rstDst As Recordset 

Say "Starting Computation of Min Costs..." 
Do 

' Find the next lowest cost vertex: 

sSQL = "SELECT Min (COST) FROM PrimerPair WHERE DONE = No" 
Set rstMin = CurrentDb . OpenRecordset (sSQL) 

If IsNull (rstMin. Fields (0) ) Then Exit Do ' Exit here! 

0 IngMinSrcCost = rstMin. Fields (0) 

"7 SSQL = "SELECT TOP 1 PPCID, DONE FROM PrimerPair WHERE DONE 

= No AND COST = " & IngMmSrcCost 

•M Set rstSrc = CurrentDb .OpenRecordset (sSQL) 

:5! IngSrcID = rstSrc . Fields ( 0 ) 

g ' Traverse all edges from that vertex: 

M If ingSrcID Mod 20 = 0 Then 

Say "Traversing all edges from " & IngSrcID 

y5 End If 

sSQL = "SELECT * FROM Edge WHERE Src = " & IngSrcID 

Set rstEdge = CurrentDb . OpenRecordset (sSQL) 

Do Until rstEdge. EOF 

' Edge cost : 

IngEdgeCost = rstEdge! COST 
' Destination: 

sSQL = "SELECT * FROM PrimerPair WHERE PPCID = " & 



m 



rstEdge !Dst 



Set rstDst = CurrentDb. OpenRecordset (sSQL) 

• See if the destination has a better path to it: 

If IngMinSrcCost + IngEdgeCost < rstDst! COST Then 

With rstDst 

.Edit 

!COST = IngMinSrcCost + IngEdgeCost 
!PRED = IngSrcID 
. Update 
End With 
End If 
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rstEdge . MoveNext 

Loop 

' Now that vertex is DONE: 
rstSrc .Edit 
rstSrclDONE = True 
rstSrc .Update 

Loop 

Set rstMin = Nothing 
Set rstSrc = Nothing 
Set rstEdge = Nothing 
Set rstDst = Nothing 
Say "Ready" 



End Sub 



5 6. FIND BEST PATH 

Jl Option Compare Database 
gi Option Explicit 

o 

|i, ' Flag PrimerPairs as SELECTED, starting at the END and working backwards, 
until there is no Predecessor: 

» Assume there is a record with REND = 10 0 000 0 000 

S'J Piiblic Sub FindBestPathO 

Dim sSQL As String 

Dim rstPP As Recordset 

Dim vntPredID As Variant 

Say "Selecting Optimal Primer Pairs..." 

sSQL = "SELECT * FROM PrimerPair WHERE REND = 1000000000" 
Do 

.Set rstPP = CurrentDb.OpenRecordset (sSQL) 
With rstPP 
. Edit 

! SELECTED = True 
.Update 
End With 

vntPredID = rstPPlPRED 

sSQL = "SELECT * FROM PrimerPair WHERE PPCID = " & vntPredID 
Loop Until IsNull (vntPredID) 
Set rstPP = Nothing 
Say "Ready" 
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End_Sub 

' New Get Cost function (starting with chr 20) : 

Private Function GetCost (ByVal IngSrcREND As Long, ByVal IngDstFPOS 
As Long) As Long 

If IngSrcREND < IngDstFPOS Then 

GetCost = (IngDstFPOS - IngSrcREND) * mcnGapPenaltyPerBase 

' Gap cost 

Else 

GetCost = IngSrcREND - IngDstFPOS 

• Overlap cost 

End If 

GetCost = GetCost + racnPenaltyPerAmplicon 

• Amplicon cost 
Q End Function 




25111369,1 
(10107967/101 IN-l) 



38 



U.S. EXPRESS MAIL #EU186312650US 



